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Abstract. The excess of continuum 7 -ray emission from the 
Galaxy above 1 GeV is an unsolved puzzle. It may indicate 
that the interstellar nucleon or electron spectra are harder than 
local direct measurements, as could be the case if a local source 
of cosmic rays were to dominate the nearby flux. It is however 
difficult to distinguish between the two cases. Cosmic-ray sec¬ 
ondary antiprotons provide a way to resolve this issue. 

We have made a calculation of the cosmic-ray secondary 
antiproton spectrum in our model, which computes self- 
consistently propagation of primary and secondary nucleons, 
and electrons. Fragmentation and energy losses are computed 
using realistic distributions for the interstellar gas and radia¬ 
tion fields, and diffusive reacceleration is also incorporated. 
Our study shows that accurate measurements of the antipro¬ 
ton flux, especially at high energies, could provide a diagnostic 
of the interstellar nucleon spectrum allowing us to test the hard 
nucleon spectrum hypothesis. Present antiproton data above 3 
GeV indicate that it can already be excluded at the few a level. 
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1. Introduction 

The spectrum of Galactic 7 -rays as measured by EGRET shows 
enhanced emission above 1 GeV in comparison with calcu¬ 
lations based on locally measured proton and electron spec- 
tra assuming the same spectral shape over the whole Galaxy 


(Hunter et al. 1997 Gralewicz et al. 1997; VIori 1997; SM97 


MS98a ). The 7 -ray observations therefore indicate that their 
spectra on the large scale in the Galaxy could be different. 
Harder cosmic-ray (CR) spectra could provide better agree¬ 
ment, but the 7 -ray data alone cannot discriminate between 
the 7 r°-decay and inverse Compton explanations (MS98b). Al¬ 
though the hard electron spectrum hypothesis seems to be more 
likely due to the probably clumpy distribution of electrons at 


high energies (e.g., Pohl & Esposito 19981, the hard nucleon 


spectrum cannot be ruled out. Explicitly, we consider the case 
that the local nucleon spectrum is not representative of the re- 
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gions within a few kpc of the sun, as could occur if a nearby 
source of cosmic rays dominates the observed fluxes. 

An important clue may be provided by secondary antipro¬ 
tons in Galactic CR produced in collisions of CR particles with 
interstellar matteiQ. These are an important diagnostic for mod¬ 
els of CR propagation and provide information complementary 
to that provided by secondary nuclei such as Be, B, and heavier 
nuclei. However, unlike secondary nuclei, antiprotons reflect 
primarily the propagation history of the protons, the main CR 
component. The observed intensities depend on the spectrum 
of CRs, their composition, details of the nuclear cross sections, 
and propagation in the Galaxy. Because they are secondary, an¬ 
tiprotons reflect the large-scale nucleon spectrum independent 
of local irregularities in the primaries. 

Previous calculations of secondary p’s have been 
made on the basis of the leaky box model (e.g., 


[Gaisser & Schaefer 1992 ; Simon & Heinbach 1996 ) and 


the locally observed nucleon spectrum. Recently several ex¬ 
periments have provided improved data on both the p/p ratio 
and the p spectrum itself (Hof et al. 19% Vlitchell et al. 1996 


Boezio et al. 1997 [Moiseev et al. 1997), and the latest calcu¬ 


lations by Simon et al. (1998) indicate good agreement with 
the data. 

We have developed a propagation code which aims to re¬ 
produce self-consistently observational data of many kinds re¬ 
lated to CR origin and propagation: direct measurements of 
nuclei, electrons and positrons, 7 -rays, and synchrotron radi¬ 
ation. These data provide many independent constraints on any 
model and our approach is able to take advantage of this since it 
must be consistent with all types of observation ([Strong 1996; 


[SM97 ; [MS98a ). In this paper we present results on the evalua 
tion of the p spectrum and p/p ratio in a model including diffu¬ 
sion and reacceleration and different nucleon injection spectra. 
Our aim is to show that p’s provide a critical test of the alterna¬ 
tive explanations of the GeV 7 -ray e xcess. O ther secondaries, 
such as positrons, also prov ide a tes t ( MS98b ), but are more af¬ 
fected by energy losses. In [MS98a we considered the positron 
fraction as evidence favouring a hard nucleon spectrum, but 


Secondary origin of CR antiprotons is basically accepted, though 
some other exotic contributors such as, e.g., neutralino annihilation 


(Bottino et al. 1991) are also discussed. 
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the spectrum considered was not as hard as required to repro¬ 
duce the 7 -ray data, and also absolute positron fluxes were not 


available at that time. In MS98b we show that positron results 


indeed confirm the conclusion of the present paper. 


2. Description of the models 

The models are three dimensional with cylindrical symmetry in 
the Galaxy, and the basic coordinates are {R, z,p), where R is 
Galactocentric radius, z is the distance from the Galactic plane, 
and p is the total particle momentum. The propagation region is 
bounded by (i?h, ±Zh) beyond which free escape is assumed. 
We take i?h = 30 kpc, Zh = 4 kpc since this is consistent with 
our B/C and ^°Be/®Be study (Strong & Moskalenko 1998a; 


SM98b). For a given Zh the diffusion coefficient as a func¬ 


tion of momentum is determined by B/C for the case of no 
reacceleration; if reacceleration is assumed then the reacceler¬ 
ation strength (related to the Alfve n speed, va) is constr ained 
by the energy-dependence of B/C ( Seo & Ptuskin 199^ . The 
spatial diffusion coefficient for the case of no reacceleration is 
taken as = l3Do{p/poY^ below rigidity po, PDo{p/poy'^ 
above rigidity po- The spatial diffusion coefficient with reac¬ 
celeration is Dxx = PDo{p/poY with ^ for all rigidities, 
and the momentum-space diffusion coefficient is related 


to Dxx ( Berezinskii et al. 1990| ; [Seo & Ptuskin 1994 ). The in¬ 
jection spectrum of nucleons is assumed to be a power law in 
momentum. The values used are Dq = 3.5 x 10^® cm^ s“^, 
Po = 5 GV, (5i = —0.60, and 82 = +0.60 for nonreaccel¬ 
eration models, and Dxx = 6 x 10^® cm^ s“^ at 3 GV and 
UA = 20 km s“^ for reacceleration models. 

The interstellar hydrogen distribution uses HI and CO 
surveys and information on the ionized component; the He¬ 
lium fraction of the gas is taken as 0.11 by number. En¬ 
ergy losses for electrons and nucleons are included (SM98b). 
The distribution of CR sources is chosen to reproduce the 
CR distribution determined by analysis of EGRET 7 -ray data 
(Strong & Mattox 1996). The secondary nucleon source func¬ 
tions are computed from the propagated primary distribution 
and the gas distribution. The 7 -ray emission from 7 r°-decay, 
inverse Compton and bremsstrahlung are computed explicitly 
in 3D from the propagated nucleon and electron spectra. 

The calculated B/C ratio is shown in Eig.|^ together with 
recent data, and the agreement indicates that our propaga¬ 
tion models are satisfactory. Our preliminary results were pre¬ 


sented in SM97 and full results for protons, Helium, positrons, 


and electrons in MS98a Evaluation of the B/C and ^°Be/®Be 
ratios, evaluation of diffusion/convection and reacceleration 
models, and limits on the halo size, as well as full de¬ 
tails of the methods are summarized in SM98b. More de¬ 


tails are available on the WWW (http://www.gamma.mpe- 
garching. mpg. de/~ aw s/aws. html). 


2.1. Antiproton cross sections 

We have used a ‘standard’ formalism to calculate p pro¬ 
duction and absorption in the interstellar medium. Antipro- 
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Fig. 1 . B/C ratio for the models with (solid) and without reacceler- 
ation (dashed). <1> = 500 MV. Data: vertica l bars: HEAO-3. Voyag er 
(Webber etal. 1996), filled circles: Ulysses (DuVemois et al. 1996 


ton production in pp-collisions has been calculated using 
the Tan & Ng (1983a) parametrization of the invariant p- 
production cross section. The total pp inelastic cross sec¬ 
tion has been calculated using a fit by Tan & Ng (1983b). 
The cross section for p production in proton-nucleus and 
nucleus-nucleus interactions has been obtained (following 
[Gaisser & Schaefer 1992 ) by scaling the pp invariant cross sec¬ 
tion with a factor = (^icrp“' + Atcr™®')/ 2 o-pp®'. where 

Apt are the atomic numbers of the incident and target nuclei. 
Eor the cross sections cr™^ and we adapted parametriza- 
tions by Tan & Ng (1983b) and Letaw et al. (1983), cor¬ 
respondingly. The p absorption cross section on an arbitrary 
nuclear target has been scaled by using the measured 
pD'^C cross section (Denisov et al. 1973; [Carroll et al. 1979 


iMakamura et al. 19^, fCuzichev et al. 1994). 

Simulations of the p production with the Monte 


Carlo model DTUNUC ( [Simon et al. 1998| ), which appear to 
be more accurate than simple scaling, have shown that He nu¬ 
clei contribute about 18% to the total p yield and their contri¬ 


bution remains a constant above the kinetic energy Tp 


500 


MeV. Heavier nuclei contribute at about the 3% level. There¬ 
fore, even if our simple scaling lowers the p yield on nuclei by 
a factor of 2 (which is unlikely at Tp ^ 500 MeV), then the to¬ 
tal yield is not underestimated by more than 10 %. In fact, other 
uncertainties dominate the secondary production, for example 
the form of the interstellar nucleon spectrum. 

Another simplification is that p’s surviving after an inelastic 
collision are totally ignored. However, calculations made with 
only the annihilation cross section show that the difference is 
small and the effect can be neglected. 


3. 7 -rays 

Eig. ^ (left) shows as an example the 7 -ray spectrum of the in¬ 
ner Galaxy for a (‘normal’) model which matches the directly 
observed electron and nucleon spectra (the latter is shown in 
Eig. Qleft). The fit to the EGRET spectra is satisfactory from 
30 to 500 MeV and the deficit above 1 GeV is evident, as dis- 
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Fig.2. Left panel: Gamma-ray spectrum of inner Galaxy (330°<i<30°, —5°<6<5°) as measured by EGRET ( Strong & Mattox 1996| ) 
compared to model with ‘normal’ nucleon and electron spectra. Also shown are the contributions of individual components: bremsstrahlung, 
inverse Compton, and 7 r°-decay. Right panel: The same compared to the model with the hard nucleon spectrum (no reacceleration). 



Fig. 3. The EGRET 7 -ray spectrum of inner Galaxy compared to 
model with the hard electron spectrum (no reacceleration). Individual 
components are the same as in Eig. ^ 


cussed in the Introduction. Simple rescaling of either electron 
or nucleon spectra does not allow the agreement to be signfi- 
cantly improved. 

A model with a hard nucleon injection spectrum (no reac¬ 
celeration, injection index 1.7) is shown in Fig. I (right). 
The corresponding propagated interstellar proton spectrum is 
shown in Fig. 0 Fig. I shows a model with a hard electron in¬ 
jection spectrum (no reacceleration, injection index 2.0). Both 
models reproduce approximately the observed spect rum, an d 
latitude and longitude profiles, almost equally well ( MS98b ), 
and hence it is difficult to discriminate between them. 


The same nucleons which contribute to the GeV 7 -ray 
emission through the decay of 7r°-mesons also produce sec¬ 
ondary p’s (on the same interstellar matter). The harder nucleon 
spectrum hypothesis, therefore, can be tested with reliable mea¬ 
surements of CR p’s. Above Tp ~ few 10 GeV the mean energy 
of parent protons is about 10 times larger than the kinetic en¬ 
ergy of produced p’s, and roughly the same holds for 7 -rays, so 


10 GeV p’s and 7 ’s both are produced by : 
Thus, the test is well tuned. 


GOO GeV nucleons. 


4. Antiprotons 

First we consider the ‘normal’ case, with nucleon injection 
spectra which after propagation and modulation match those 
locally observed (Fig. ^ left). Our calculations of the interstel¬ 
lar p spectra and p/p ratio for these spectra are shown in Fig. ^ 
The computed p spectrum is divided by the same interstellar 
proton spectrum, and the ratio is modulated to 750 MV. The 
corresponding ratios are shown on the right panel. We have 
performed the same calculations for models with and without 
reacceleration and the results differ only in details. As seen, our 
result agrees well with the calculations of Simon et al. (1998), 
showing that our treatment of the production cross-sections is 


adequate as discussed in Section 2.1 


We now turn to the case which matches the 7 -ray data at the 
cost of a much harder proton spectrum than observed (Fig. ^ 
right). The dashed lines in Fig. ^ (right) show the p/p ratio for 
the hard proton spectrum (with and without reacceleration); the 
ratio is still consistent with the data at low energies but rapidly 
increases toward higher energies and becomes ~4 times higher 
at 10 GeV. Up to 3 GeV it does not confict with the data with 
their very large error bars. It is however larger than the point at 


3.7-19 GeV ( pTof et al. 1996| ) by about 5cr. Clearly we cannot 

, but it does 


conclude definitively on the basis of this one ] 
indicate the sensitivity of this test. In view of the sharply rising 
ratio in the hard-spectrum scenario it seems unlikely that the 
data could be fitted in this case even with some re-scaling due 
to propagation uncertainties. It is interesting to note that the 
local p/p ratio seems to depend only slightly on the details of 
the propagation. 

^ We do not consider here the older p measurement of Golden et al. 
(1984) b ecause the flight of the early instrument in 1979 was repeated 
in 1991 ( Hof et al. 1996 ) with significantly improved instrument and 
analysis techniques. Thus the latter data are more reliable and the rel¬ 
evance of this measurement to the earlier one is indicated in Hof etal. 
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Fig. 4. Left panel: Interstellar nucleon and antiproton spectra as calculated in nonreacceleration models (thin lines) and models with reaccel¬ 
eration (thick lines). Spectra consi stent with the loc al one are shown by the solid lines, hard spectra are shown by the dashed lines. The local 
spectrum as measured by IMAX (VIenn et al. 1997) is shown by dots. Right panel: p/p ratio for different ambient proton spectra. Lines are 
coded as on the left. The ratio is modulated with 4? = 750 MV. Calculations of Simon et al. (1998) are shown by the dotted lines. Data from: 
■ Boezio et al. (1997), O Bogomolov et al. (1987,1990), A Hof et al. (1996), □ Mitchell et al. (1996), 0 Moiseev et al. (1997). 


Our main conclusion is that antiprotons provide a sensi¬ 
tive test of the interstellar nucleon spectra and hypotheses for 
the origin of diffuse Galactic 7 -rays. On the basis of the p/p 
data point above 3 GeV we seem already to be able to ex¬ 
clude the hypothesis that the local CR nucleon spectrum differs 
significantly from the Galactic average (by implication adding 
support to the ‘hard electron’ alternative), but confirmation of 
this conclusion must await more accurate data at high ener¬ 
gies. In this respect we note that the p/p ratio from Hof et al. 
(1996) is currently being refined and absolute p fluxes will be 
calculated (Hof 1998, private communica tion). Additionally , 
a re-flight of the CAPRICE instrument ( Boezio et al. 1997 ) 
took place in spring 1998, and several other balloon instru- 
ments could be ada pted for anti p roton measurements (H EAT: 
Barwick et al. 1997), ISOMAX: Streitmatter et al. 1993). On 


longer timescale several satellite e xperiments are plan ned or 
under constructio n (e.g., PAMELA: Adriani et al. 1995 ; AMS: 
Ahlen et al. 1994 ). These new experiments should allow us to 
set stricter limits on the nucleon spectra including less extreme 
cases than considered here, and to constrain better the interpre¬ 
tation of 7 -rays. 
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